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SATURATION SCALE FROM THE BALITSKY KOVCHEGOV 

EQUATION 

KRZYSZTOF GOLEC-BIERNAT 
H. Niewodniczanski Institute of Nuclear Physics 
Radzikowskiego 152, 31-342 Krakow, Poland 
E-mail: golec@ifj.edu.pl 

Recently Munier and Peschanski presented an analysis of the Balitsky-Kovchegov 
(BK) equation concerning the extraction of the saturation scale, using a simpler 
equation. We numerically analyze the full BK equation confirming the universality 
of their analysis. 

1 Introduction 

The Bahtsky-Kovchegov equation lll2j describes saturation 3 in gamma - nu- 
cleus deep inelastic scattering, when the structure function F2 ~ log(l/a;B) for 
the Bjorken variable a:^ — *■ 0. To be more precise, this behaviour appears in the 
region where the photon virtuality < Q^{xb)- Here is the saturation scale 
which rises with decreasing xb- Thus the famous BFKL growth 4 in the region 
> QI{xb), where F2 ~ Xg'^, is tamed due to saturation effects generated by 
the BK equation. The emergence of the saturation scale and its energy dependence 
is the key issue in analyses of saturation. Recently, a very elegant method of the 
extraction of the saturation scale was developed in 0, based on a simplified ver- 
sion of the BK equation. According to mathematical properties of this equation, 
the obtained results are universal and concern also the full BK equation. In this 
presentation we show numerical results which support this claim. 

2 Munier Peschanski analysis 

In the space of transverse gluon momentum fc, the BK equation reads 

dycj) ^asx{-dL)(l> ~as(f>'^ , (1) 

where (j>{Y, L) is the Fourier transform of the quark dipole-nucleon scattering am- 
plitude Here x{l) = 2V'(1) — "0(7) — '0(1~7) is the BFKL characteristic function 
13], rapidity Y — \n{l/xB), and L = ln(fc^/A;g). The structure function F2 at small 
Xb can be computed from solutions of the BK equation . 

Eq. ^ is nonlinear, with the linear part given by the BFKL function. The 
nonlinearity is responsible for saturation, taming the BFKL growth generated by 
the linear part alone. In a remarkable simplification was done approximating 
the hnear part by the Taylor expansion xil) — Xc + Xcil - 7c) + ^Xc(7 - 7c)^, 
where 7c — 0.63 obeys the equation jcX'c = Xc- Then the linear change of variables 
from (y, L) to the new ones {t, x) was performed 

t = i (a.7c' Xc )y a; = 7c L + (7,2 x'^ - Xc) F, (2) 
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together with the redefinition: u(t,x) = 2 / (j'^Xc) 4'{Y, L) . As a result, the Fisher- 
Kolmogorov-Petrovsky-Piscunov (FKPP) equation was found, 

dtu^dlu + u{l-u). (3) 

Eq. |(2Jl admits traveling wave solutions: 

>00,x) ~ u{x — 171/3(1)), (4) 

where the function mp{t) depends on the asymptotic behaviour of the initial con- 
dition u{tQ, X oo) ~ exp{— /3a:}. For example, for (3 > 1 

m(i{t)^2t-l\nt + 0{l). (5) 

This case is particularly important since it contains initial conditions with color 
transparency, 4> ^ l/fc^ for k —> oo, corresponding to /S — l/jc- From the form (0J, 
the function has the property of geometric scaling 6 , 

c^{Y,k)^c^{k^/Ql{Y)), (6) 

with the saturation scale given by 

InQliY) ^^Y-^lnY~^.^^^+0(l/Y), (7) 
7c 7c'V "-^"VF 

where the third term was computed in The first three terms are the only 
universal terms since the shift Y — > Y + Yq, which amounts to the change of an 
initial condition, affects the contribution 0{1/Y). 



3 Numerical analysis 

The purpose of the numerical analysis is to confirm that the results ^ and ^ are 
also valid for the full BK equation. To this end we numerically solved eq. ^ with 
the help the Chebyshev polynomial method. 

In Fig. 1 we show the solutions 4>(Y, k) as functions of transverse gluon momen- 
tum k for twenty different values of rapidity from Iq = to F = 40. The two plots 
show the transition from the power-like behaviour (f) ~ for large k (the r.h.s. 

plot in the log-log scale) to the logarithmic behaviour ^ In fc for fc ^ (the same 
plot in the log-lin scale on the l.h.s). This is the illustration of the transition to satu- 
ration due to the nonlinear term in eq. , when the power-like growth of the gluon 
distribution for decreasing k is eventually tamed and logarithmic behaviour sets in. 
The initial condition (the dashed lines) contains color transparency and saturates 
to a constant value for fc ^ 0. The latter behaviour is immediately changed by the 
nonlinear evolution of the BK equation to the logarithmic dependence on k. The 
FKKP equation, however, conserves the saturation value of the initial condition 
producing the traveling wave moving to the right, see 0]. Thus the truncation of 
the BFKL kernel in the Munier-Peschanski analysis influences only the saturation 
region while for large k the two solutions are identical. 
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Figure 1. The solution (j)(Y,k) of eq. (1) for twenty values of rapidity from Yq = to y = 40 in 
the order indicated by the arrow. These are the same plots in the log-lin scale (a), and the log-log 
scale (b). They illustrate the transition to saturation with decreasing k, when the power-like 
behaviour (b) changes into the logarithmic one (a). The dashed lines is the initial condition with 
color transparency, <f> ^ 1/fe^ for large k. 



How to extract the saturation scale from numerical solutions ? The logarithmic 
pattern shown in the Fig. 1(a) suggests that after an initial period of the evolution, 
when details of the initial condition are washed out, the following behaviour sets in 

mk)=\n{Qs{Y)/k). (8) 

This could be checked numerically if the saturation scale computed from ©, 
Qs{Y) — fcexp{0(y, /c)}, is independent of k. For this purpose we compute Qs{Y) 
for four different values of logj^Q k — —2, —1,0,1. In this range of k excellent agree- 
ment was found, which is shown in Fig. 2(a) where the four saturation scales coin- 
cide. All momentum scales are in units of some arbitrary scale ko, thus the absolute 
normalization of the saturation scale is not determined. In order to compare the 
y-dcpendencc of our saturation scale with the dependence given by eq. iQ, we 
compute the derivative dluQ^/dY, which is independent of the absolute normal- 
ization of Qs- The result is shown in Fig. 2(b). The solid lines correspond to the 
derivative computed from eq. Q with one, two or three first terms in the asymp- 
totic expansion. The lower dashed line is given by differentiation of the saturation 
scale defined by eq. ||HJ). We see that such a numerically computed Qs agrees with 
the expansion Q with the first two terms. The third term is still not visible at the 
maximal rapidity achieved in the numerical analysis. It remains to be answered 
whether this is a numerical effect or the agreement is achieved at higher values of 
rapidity. 

There is another method of computation of the saturation scale, used in most 
of the analyses based on the linear BFKL equation, i.e. from the equation 
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Figure 2. In figure (a) the saturation scale computed from is shown for four values of k given 
in the text with no visible difference between them. In figure (b) the derivative d\nQg/dY is 
shown. The solid lines correspond to the derivative of the analytic formula J7J with (as indicated) 
one, two or three terms. The lower dashed line is computed from the saturation scale defined by 
eq. jSJ, while the upper dashed line is computed from formula j5J. 



This method is more appropriate for the comparison with the Munier-Peschanski 
analysis since it uses the solution in the wave front region. In Fig. 2(b) the upper 
dashed curve corresponds to the saturation scale computed from eq. 0, which 
agrees very well with the three term formula 0. 

In conclusion, the analysis of the geometric scaling and saturation scale in the 
FKKP equation is general and applies to the full nonlinear BK equation. 
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